States of fermionic atoms in an optical superlattice across a Feshbach resonance 



T. Goodman and L.-M. Duan 

FOCUS center and MCTP, Department of Physics, University of Michigan, Ann Arbor, MI 48109 

We investigate states of fermionic atoms across a broad Feshbach resonance in an optical super- 
lattice which allows interaction only among a small number of lattice sites. The states are in general 
described by superpositions of atomic resonating valence bonds and dressed molecules. As one scans 
the magnetic field, level crossing is found between states with different symmetry properties, which 
may correspond to a quantum phase transition in the many-body case. 
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I. INTRODUCTION 

The past several years have seen many new develop- 
ments in ultracold atomic physics Q, ■ Two experimen- 
tal control techniques have played a critical role in these 
advances: the use of optical lattices to produce diverse 
interaction configurations , and Feshbach resonance to 
control the strength of the interactions between atoms 
Q . This has motivated significant interest in combining 
these two techniques MM H H EJ 

Close to a broad Feshbach resonance, one then has 
strongly interacting fermions in a lattice, with the effec- 
tive interaction energy easily larger than the band gap. 
For such a strongly correlated system, in general it is very 
hard to understand its physics. As the first step, one may 
try to understand the single-site physics, which contains 
the dominant interactions for the atoms in a deep optical 
lattice. The two-body physics on a single site has been 
solved exactly in [j| under typical experimental configu- 
rations, and it is found there that many bands of the lat- 
tice are significantly populated, in agreement with the ex- 
perimental observation 0. In the next step, it is natural 
to consider interactions over a few lattice sites. We have 
shown in [f| that for the case of a multiple-site lattice, one 
needs to take into account the direct collision interactions 
between the neighboring sites, as the magnitude of this 
interaction is larger than the atom tunneling rate for a 
typical wide Feshbach resonance such as with 40 K or 6 Li 
atoms. An effective lattice interaction Hamiltonian has 
been derived there, which takes into account both the 
multi-band populations and the direct off-site and on- 
site collisions. In this effective Hamiltonian, whenever 
two atoms come to the same lattice site, they form a 
dressed molecule state (a single-site Cooper pair), which 
corresponds to an exact eigenstate constructed from the 
single-site physics. The effective Hamiltonian then de- 
scribes the interaction between these dressed molecules 
and the atoms over different lattice sites. The explicit 
form of the Hamiltonian is as follows 6] 

Heff = 5>(BHjd 4 + ]T t d P4djP 
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where and represent creation operators for 
fermionic atoms and bosonic dressed molecules, respec- 
tively; <T =f, | labels two internal spin states; and i labels 
the lattice sites (with j £ N(i) labeling the sites adjacent 
to i). P represents a projection of the state at every lat- 
tice site i onto the four-dimensional subspace with basis 

||0) i , a| T |0) - , atj \0). , d\ The parameters g, t d , t a , 

and tda depend on the physical properties of the atomic 
system as well as the multi-band properties of the optical 
lattice (see the explicit expressions in 6]). The param- 
eter A corresponds to an eigenenergy of the two-body 
physics from a single site, and can be tuned over a wide 
range of values by varying the applied magnetic field B 

In this paper, we make use of the above effective Hamil- 
tonian H e ff to study the physics of this strongly corre- 
lated system with interactions among a few lattice sites. 
In particular, we focus on investigation of the states 
of a single plaquette, which is a basic unit of the two- 
dimensional square lattice. This study has two purposes. 
First, understanding the states of atoms at a single pla- 
quette is a necessary step towards the challenging goal 
of understanding the physics of this strongly interacting 
gas in a quasi- two-dimensional optical lattice. It is shown 
m that the effective Hamiltonian H e ff reduces to the 
well-known t- J model [Tj and the XXZ model (Tj] in cer- 
tain limits of the parameter values, so the physics associ- 
ated with H e ff should be rich and parameter-dependent. 
In this work, we would like to understand the influence 
of the parameters on a few body physics, and that under- 
standing will provide an intuition for taking appropriate 
approximations towards the many-body physics. In par- 
ticular, from the few-body physics, we construct states 
which will provide the basic entries for an effective many- 
body theory through the contractor renormalization pro- 
cedure (a real-space renormalizaton group method for 
high dimensions) |lall4| . We will see that even for a sin- 
gle plaquette, the behavior of the states has been pretty 
rich. The eigenstates involve resonating valence bonds 
and superposition of dressed molecules, and are highly 
entangled over different lattice sites. With variation of 
the parameters in the Hamiltonian -ff e //, there are sev- 
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FIG. 1: Superlattice potential vs. x for V(x) — 
_(V lS m 2 (If) + V 2 sin 2 



eral level crossings for the lowest eigenstate with change 
of the state symmetry properties, which may correspond 
to a quantum phase transition for larger systems. 

Second, the study of the physics of a single plaque- 
tte is also of practical relevance. For atoms in an opti- 
cal superlattice, the physics can be dominated by the 
interactions within single plaquettes. A simple opti- 
cal superlattice can be formed by adding two stand- 
ing wave laser beams with commensurate wave vectors 
0,0. The potential, say, along the x direction, has 
the form V(x) = — [Visin 2 {k\x) + V^sin 2 (fe^)] ■ If we 
choose the wave vector &2 = 2fci = 2ir/L and require 
< V\ < 4V2, then we have potential barriers of two 
different heights (see FIG. The minima occur at 
x = nL ± xq (for integer n) , where 

- - 5"-' @) < 2 » 

The lower and higher potential barriers Vi ow and Vhigh 
are given respectively by 

V hiah = (3) 

The barrier Vhigh can be significantly larger than V\ ow if 
we choose V\ close to AV2, and such a high barrier turns 
off the interactions except for the ones between the sites 
separated by Vi ow . If we apply this optical superlattice 
potential along both the x and y directions and a deep 
lattice potential along the z direction, we then have inter- 
actions dominantly within the single plaquettes in the x-y 
plane. With strongly interacting atoms in this optical su- 
perlattice potential, one can test the predictions from the 
effective Hamiltonian H e 1 1 , and detect the exotic entan- 
gled states emerging from the ground state configurations 
of H eff . 

We should also point out that the effective Hamiltonian 
H e ff includes the Hubbard model as a particular case. 



The Hubbard model is given by the Hamiltonian 0, 0] 
H Hu b = -t ^ ( a L a j° + H - c -) + U *^2 ni T"-ii> ( 4 ) 

<i,j>,a i 

where — a^a^. Specifically, H e ff can be written 
in the form of Huub if we substitute d\ with a\^a\^ and 
make a particular choice of the parameters in H e f f with 
t a = —t, tda = t, g = t, td = 0, and A = U. So, one can 
see that H e ff extends the well-known Hubbard model 
Huub m a nontrivial way. Note that for strongly interact- 
ing atoms near a broad Feshbach resonance, the parame- 
ters g and td a are significantly different from the atomic 
tunneling rate t due to the multi-band populations and 
the direct neighboring collisions. From the expressions of 
these parameters in Ref. |q, we estimate that typically 
\td\ « \ta\ « \tda\ ~ |fli- This is because t a corre- 
sponds to atomic tunneling in the single lowest band, 
whereas td a and g correspond to interactions involving 
the dressed molecule states (which are superpositions of 
states in multiple upper bands). For the numerical cal- 
culations in this work, we typically take td ~ 0, — t a ~ 
0.1|<?|-0.3|g|, and tda ~ |ff|-2|<?|. (Note from the form of 
H e ff that the sign of g is essentially irrelevant, as it can 
be incorporated into the definition of d' .) The parameter 
A is sensitive to the external magnetic field, and can be 
scanned from the value much smaller than — \g\ to the 
value much larger than \g\. 

The atomic states within each plaquette critically de- 
pend on the atom number and the spin configuration in 
that plaquette. In the following, we will consider all the 
different nontrivial cases with different numbers of spin 
t and J, atoms occupying the four-site plaquette. 

II. FOUR ATOMS PER PLAQUETTE: TWO |, 
TWO I 

Over most of the typical range of the parameter values, 
the plaquette occupied by two f and two J, atoms has 
two distinct types of ground states, with a level crossing 
occurring at some critical value of A. These two types of 
states can be distinguished by how they transform under 
a 90° rotation in the plane of the plaquette. Under such 
a rotation, the ground state wavefunction for A less than 
(greater than) the critical value is multiplied by a factor 
of +1 (—1). Thus, we say that the phase on the negative 
side of the transition has s-wave symmetry, and that on 
the positive side has d-wave symmetry. 

The ground states of each of these two types change 
smoothly with changes in the parameter A. Thus, we 
can identify particular energy eigenstates as the "s-wave 
state" and the "d-wave state" over the full range of A, 
even as the exact form of the eigenstate changes. (Note 
that these are not the only eigenstates with s-wave and d- 
wave symmetry - here we use these terms to refer solely 
to those states which are the ground states when the 
system is in the corresponding parameter regions.) The 
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FIG. 2: Energy vs. A for a plaquette occupied by two \ and 
two J, atoms. Other parameters are td. a = l-5|p|, t a — — 0.2|g|, 
td — (a): Eigenenergies of the s-wave (o) and d-wave (x) 
states, (b): Energy difference (gap) between ground state 
and first excited state. The gap vanishes at the level crossing 
point. Because the eigenenergies vary smoothly with A, the 
curve is smooth except at the level crossing points for the 
ground state (where the gap is zero) and for the first excited 
state (indicated by arrows) 



energies of the s-wave and d-wave states can be easily cal- 
culated through exact diagonalization, and they are plot- 
ted in FIG. Eta), which illustrates the crossover between 
them. (For this hgure, we scan A and set the other pa- 
rameters of H e ff to their typical values with td a = l-5|si 
and t a = — 0.2|g|.) The energy gap between the ground 
state and 1st excited state is shown in FIG.EJb)). 

To understand the properties of the ground state, it is 
important to have its explicit expression. Although one 
can calculate this explicit expression through numerical 
exact diagonalization, the state is in general a superposi- 
tion of many basis- vectors (36 vectors in this case), with 
all the superposition coefficients varying with A. It is 
troublesome to understand the state's properties from 
this lengthy expression. To overcome this problem, we 
describe the s-wave and d-wave states more compactly, 



in a way that illustrates their rotational symmetry, by 
means of a pictorial representation which we define here. 

1 2 

The four sites of a plaquette we label as: 3 4. We place 
various pictures on these sites corresponding to creation 
operators applied at those sites. The whole picture repre- 
sents the product of these operators, applied to the vac- 
uum state |0). For instance, I placed on two sites (either 
horizontally, vertically, or diagonally) represents a nor- 
malized singlet between those two sites. So, if the sites 

are labeled i and j, this represents ( a iT a Ji — a ll a Jt) ' 
(Note that the order of i and j does not matter, as the 
anti-commutation of a\ and Oj makes the singlet sym- 
metric under exchange of i and j.) © represents a 
dressed molecule. (I.e., if located at site i , this pic- 
ture corresponds to a[.) □ represents an unoccupied site. 
The creation operators that make up a singlet are always 
grouped together; other than that, the order of the oper- 
ators is irrelevant, as the singlets and dressed molecules 

T © 

commute. As an example, the picture 1 o represents 

^ (aJ T 4 ; - ai^y 4 |0>- 

The full Hilbert space for two | and two { atoms on 
a plaquette is 36-dimensional. However, the s-wave state 
(over the full range of A) can be conveniently expressed 
as a vector in a 4-dimensional subspace of the full space, 
with basis vectors: 
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The d-wave state (as A varies) can be written as a 
vector in a 3-dimensional subspace with basis vectors: 
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Note that |4) s and |3) d are written here as sums of 
non- orthogonal terms; however, this form makes their 
rotational symmetry readily apparent. The states |4) 
and |3) d are the resonating valence bond (RVB) states 
for atoms on a sing le plaquette [H with s and d 
wave symmetries, respectively. For a larger lattice, the 
RVB states are in general superpositions of many differ- 
ent spin-singlet distribution patterns |18| . 
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FIG. 3: The ground-state configuration vs. the detuning A 
for a plaquette occupied by two f and two J. atoms, (tda = 
l-5\g\, t a — — 0.2 1 p | , td = 0) (a) Components of the s-wave 
state (si: o, S2: x, S3: y, S4 : +)• (b) Components of the 
d-wave state (di: o, efe: x, ^3: +). The marked data points 
were computed from the full Hamiltonian H e //> whereas the 
solid lines were computed from the projected Hamiltonians 
H s and Hd, respectively. 



Thus the s-wave and d-wave ground states, respec- 
tively, are: 



si 
di 



d2 \2} d 



S3 |3} s - 
- da |3> d 5 



(5) 
(6) 



They are superpositions of many different distribution 
patterns of the dressed molecules and the atomic valence 
bonds (spin singlets). The values of the superposition 
coefficients are shown in FIG.|3]as a function of the ratio 
A/|<7|. Note that in the limiting case A/|</| >> 1, the 
effective Hamiltonian H e ff reduces to the t-J model 0. 
Indeed, one can see from FIG. |3] that the state \ip) d in 
that limit tends to the ground state |3) d of the t-J model 
on a plaquette [Til ]. 

Projected onto these subspaces, H e ff (with t d — 0) 
expressed in terms of the bases shown above becomes: 



FIG. 4: In the typical range of t da /\g\ and t a /\g\, the s-wave 
triplet ground state of the plaquette with 2 f and 2 1 atoms 
occurs for parameter values within the shaded region. 
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for the s-wave state, and: 
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for the d-wave state. The lowest energy eigenstates of 
these two Hamiltonians are the s-wave and d-wave states 
(respectively) of the full Hamiltonian H e ff. (See solid 
lines on FIG. El) 

For a small portion of the range of the parameter values 
tda/\g\ and to/lsi there is an additional type of ground 
state, which occurs for A between the s-wave and d-wave 
states above. For this type, the ground state also has 
s-wave rotational symmetry. However, the states f and 
I are in a triplet configuration, rather than the singlet 
occurring in the other two types of states \if}) s and \if>) d . 
The region of the parameter space for which this triplet 
phase occurs is shown in FIG. 01 The eigenenergies of the 
s-wave singlet, s-wave triplet, and d-wave singlet states 
are shown in FIG. OJa) for t da — 2\g\, t a = — 0.3|<7| (which 
is within the range where the triplet ground state occurs.) 
The gap between the ground state and first excited state 
for t da = 2\g\, t a = -0.3|g| is shown in FIG.^b). 

The s-wave triplet state can be written as a linear com- 
bination of three states: 
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FIG. 5: Energy vs. A for a plaquette occupied by two f and 
two I atoms. Other parameters are tda = 2|gr|, t a = — 0.3|g|, 
td — (a): Eigenenergies of the s-wave singlet (o), d-wave 
singlet (x), and s-wave triplet (+) states, (b): Energy differ- 
ence between ground state and first excited state. Crossovers 
in the first excited state are indicated by arrows. 

i 3 w=i7i(n + i! + ^; + ii) 

Here #=0 represents the triplet ( a iT a ii a ll a it)' 
where i is the site of the black-filled circle, and j is the 
site of the white-filled circle. Note that unlike the singlet, 
the triplet is not symmetric under exchange of i and j: 

71 ( a !Xi + °Mt) = -js { a h a l + sVIt) 

Projected onto the basis {\l) trip , \2) trip , |3) trjp |, 
He// (with td ~ 0) becomes: 

/ _ A \/2 (**,-*») -2^2ff\ 

Htrip = V^^da-ia) A (9) 

\ -2v/2ff / 

and the ground state of this Hamiltonian is the s-wave 
triplet state. It should also be noted that the s-wave 
triplet state is the first excited state of H e f f in the limit 



FIG. 6: Energy vs. A for a plaquette occupied by three j and 
one I atoms. Other parameters are tda ~ 1-5|<?|, t a = — 0.2|p| , 
td — (a): Ground state energy, (b): Energy difference 
between ground state and first excited state. 



of large positive A. 



III. FOUR ATOMS PER PLAQUETTE: THREE 

T, ONE I 

The plaquette occupied by three f and one j atoms 
has only one type of ground state over the full range of 
A. The ground state has s-wave symmetry (i.e. it is un- 
changed under 90° rotations in the plane of the lattice). 
The ground state energy is plotted in FIG.UJa), and the 
energy difference between the ground state and the first 
excited state is shown in FIG. E{b). For this figure the 
other parameters were tda = 1 -5 1 (7 1 and t a = — 0.2\g\. 

The ground state can be represented compactly in 
the pictorial representation introduced above. Here we 
add an additional symbol to represent a single atom 
in the | state. Because the order of the fermionic 
creation operators matters, we use f to represent the left 
creation operator and to represent the right creation 
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FIG. 7: Components of the ground state (Ci: o, C2: x, C3: 
+) vs. A for a plaquette occupied by three f and one j atoms, 
(ida = 1 -5| <7 1 , to = — 0.2 1 (7 1 , td = 0) The marked datapoints 
were computed from the full Hamiltonian H e ff, whereas the 
solid lines were computed from the projected Hamiltonian 

© o 

operator. For instance, L |0), whereas 

© ° 

1r t =44 t 4tI°)- 

The ground state is: |^) 5 = Ci |1) 5 +C 2 |2> 5 + C* 3 |3) 5 , 
where 




The values of the coefficients C\, C2, C3 are shown in 
FIG.EI 

Projected onto the three-dimensional subspace of the 
full Hilbert space with basis vectors 1 1)^ , |2) s , |3) s , H e ff 
(for td — 0) is: 

/ _ A V2(t a -t da ) -2s/2g\ 

H s = V2(t a -t da ) A (10) 

\ -2y/2g / 

Thus, the ground state of this Hamiltonian is the ground 
state of H e ff. (See solid lines in FIG.0) 

It should also be noted that the s-wave ground state for 
3 t, 1 1 atoms per plaquette is degenerate with the triplet 
state for 2 | and 2 J, atoms described in the previous 
section. In fact, the 2 f, 2 J. triplet state is identical to 
the ground state for 3 f and 1 | atoms and for 1 f and 3 

I atoms, except that the triplet ( a iT a ]l ~^ a ll a ]T] ^ s 
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FIG. 8: Energy vs. A for a plaquette occupied by one f and 
one | atom. (t da = 1.5|g|, t a = -0.2|g|, t d = 0.) (a): Ground 
state energy (b): Energy difference between ground state and 
first excited state. The curve is smooth except at a crossover 
in the first excited state (indicated by an arrow). 



replaced with aLaL in the 3 T, 1 I case, and with a],a - 
in the 1 J, 3 J. case. 



IV. TWO ATOMS PER PLAQUETTE: ONE f, 

ONE I 

When occupied by only a single atom of each spin 
state, the plaquette has a single type of ground state 
for all values of A (for values of the other parameters 
within the typical range). This state is symmetric under 
90° rotations - i.e., it has s-wave symmetry. The ground 
state energy of this system is plotted in figure IHJa). Fig- 
ure IHb) shows the excitation gap between the ground 
state and first excited state. Both these figures assume 
typical values of t da and t a (t da = 1.5\g\, t a = -0.2\g\). 

The ground state can be expressed as a vector in a 3- 
dimensional subspace of the full Hilbert space. The basis 
vectors of this subspace (in the pictorial representation 
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FIG. 9: Components of the ground state (ci: o, C2: x, C3: 
V) vs. A for a plaquette occupied by one f and one J. atom, 
(ida = 1 -5| <7 1 , t a = — 0.2 1 (7 1 , td = 0) The marked datapoints 
were computed from the full Hamiltonian H e ff, whereas the 
solid lines were computed from the projected Hamiltonian H. 



introduced above) are: 

|i> = ! ( ® + ®+ : + 

> o o 00 o (rri) 

l 2 > = I ( 
I 3 ) = 75 ( X - 

Thus, the ground state is given by: \ip) = c\ |1) + 
c-2 1 2) + C3 |3), where the values of the coefficients Ci, C2, 
C3 are shown in FfG. El 

Projected onto this subspace, H e ff (for td — 0) ex- 
pressed in the above basis becomes: 
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(11) 



Thus, the ground state of this Hamiltonian is the ground 
state of H e ff. (See solid lines on FIG. El) 



V. THREE ATOMS PER PLAQUETTE:TWO T, 

ONE I 

The plaquette with two t atoms and one J, atom has 
three distinct types of ground states for different val- 
ues of the parameter A (with the other parameters in 
the typical range). However, over a wide range of A 
around A = the system is in the same type of ground 
state. The ground state of this type is two-fold degen- 
erate. (Hence, we will refer to this as the "degenerate 
state".) The ground state energy and the gap between 
the ground state and first excited state are shown in FIG. 
ITU1 (For the full range of A values shown in the figure, 
the system is in the degenerate state.) The degenerate 
ground states (which we call \4>) + and |?/')_) can be de- 
fined in such a way that they are eigenstates of a 90° 




FIG. 10: Energy vs. A for a plaquette occupied by two | 
and one 1 atoms. (t da = 1.5|g|, t a = — 0.2|g|, t d = 0) (a): 
Ground state energy (b): Energy difference between ground 
state and first excited state. Crossovers in the first excited 
state (at which points the curve is not smooth) are indicated 
by arrows. 



rotation in the plane of the plaquette, in which case \?p)± 
gains a factor of ±z under such a rotation. 

The state \ip} + can be expressed as a vector in a par- 
ticular six-dimensional subspace of the full Hilbert space. 
We define the basis vectors of this subspace in the picto- 
rial representation introduced above. However, because 
the order of the fermionic creation operators matters, we 
use three symbols |, ff, A Q, JJ-, V) to represent the first, 
second, and third creation operator for atoms in the | 
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(|) state. E.g., -U- 1 = a^a^a^ 
this way, the six basis vectors are: 
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and the state |^>) + is given by: 
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|V>, = A|1) + + A*|2) + +B|3) + + C|4K 



(12) 



for some complex coefficients A, B, C, and D. Note 
that under a 90° clockwise rotation \ri) , — > i |n) , for 
each n, and thus 1^) , — ► i |^) , . The state \4>)_ can be 
expressed as a vector in a six-dimensional subspace of the 
full Hilbert space with basis vectors 

|n)_ = \n)l (13) 

for n = 1, 6. _ is given by: 

\i>)_ = A*|1)_+A|2)_+.B*|3)_ + C*|4)_ 

+D*|5)_+D|6>_ = |</>>; (14) 

The complex coefficients A, f?, C, and D can be written 
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S = |B|e l7r/4 , C= |C|e l7r/4 



(15) 



where |C|, |D|, 0^, and D depend on the pa- 

rameters of H e ff. The coefficients are determined up to 
an arbitrary overall phase, which here was chosen to fix 
the phases of B and C as shown. (B and C were found to 
have the same phase.) For typical values of t a = —0.2\g\, 
tda = l-5|gi, td — 0, the values of these parameters (vs. 
A) are shown in FIG. El 

Projected onto the subspace with basis 
{|1) + ,|2> + ,|3) + ,|4) + ,|5> |6) + }, H eff (expressed 
in that basis) is (for td — 0): 
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The ground state of H + is thus 1^)+- (See solid lines 
in FIG. El) Projected onto the subspace with basis 
{|1)_,|2)_,|3)_,|4)_,|5)_,|6)_}, H eff is given by: 

H-=H* + = Hi (17) 

(Note that in this equation H + is still expressed in the 
basis in which it was defined above.) Thus the ground 
state of H_ is \tjj}* + = 

For A far to the negative side (A < —92.9 \g\ for 
tda = t a = -0.2\g\, t d = 0), the system of two 
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FIG. 11: Ground state parameters vs. A for a plaquette 
occupied by two \ and one \ atoms. (td a — 1-5|<?|, t a = 
-0.2|g|, t d = 0) fa): Amplitudes (\A\: o, \B\: x, \C\: v, 
\D\: +). (b): Phases (4>a'- °, 4>d'- x). The overall phase 
was chosen to give B = \B\e iw/4 and C = \C\e iw/4 . The 
marked datapoints on (a) and (b) were computed from the 
full Hamiltonian H e f / , whereas the solid lines were computed 
from the projected Hamiltonian H+. 



\ atoms and one J. atom on a plaquette has a non- 
degenerate d-wave ground state. This state can be ex- 
pressed as a vector in a 3-dimcnsional subspace of the 
full Hilbert space of this system, with basis vectors: 
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|3),«/«=K y.-y + y.-y.) 
= H i ! - ~ + ; i - T „ 
- 1 ; + n - ! i + -j.) 

Note that there is some ambiguity in which spins to 
group into a singlet and which to write as f, as evidenced 
by the two forms given for |3); e y t (the first of which has 
the advantage of all its terms being orthogonal, but the 
second of which makes more obvious how the Hamilto- 
nian connects it to components |l); e ^ t and |2) Je y t ). 

H e ff (for td — 0) projected onto this subspace ex- 
pressed in the basis ||l) ie/t , \2) left , |3) ie/t } is: 

/ A-tda -V2t a -g \ 
Hieft = -V2t a A -V2g (18) 
V -g -V2 5 t a J 

Thus, the ground state of this Hamiltonian is the ground 
state of H e ff in the left-most region (A < —92.9 \g\). 

For A far to the positive side (A > 97.9 \g\ for td a = 
1.5\g\, t a = — 0.2\g\, td — 0), the system has a non- 
degenerate s-wave ground state. Furthermore, in this 
state the ground state wavefunction and energy are con- 
stant for changing A. In the pictorial representation this 
ground state is given by: 

+ n + 1 i + a + i T * 

+ X + X + V. + X) 

The energy of this state is H r i g ht = 2t a . 

It should be noted that the case of 1 |, 2 J, atoms on a 
plaqucttc is equivalent to the 2 "f , 1 J, case under exchange 
of | and J, spins. The Hamiltonian H e ff is invariant un- 
der such a spin exchange, except for a change in the sign 
of g. This is equivalent to replacing d' with — dJ. Thus, 
the eigenenergies of these two cases are identical, and the 
eigenstates are identical except for a change in the sign 
of the components which include a dressed molecule. 

VI. SUMMARY AND DISCUSSION 

In the above, we have investigated the ground state 
properties of the system with different numbers of spin 
t and I atoms occupying the four-site plaquette in an 



optical superlattice. All the other cases can be reduced 
to one of the configurations considered above, or to a 
trivial case, through the particle-hole exchange. (Cases 
where all particles are in the same spin state are non- 
interacting, and thus trivial.) For instance, for five atoms 
with three spin- J. and two spin-| , one has two spin-f and 
one spin-| holes in that plaquette. So, the states are 
equivalent to those in the case with two spin-t and one 
spin-| atoms, but with exchange of the parameters tda 
and t a in the effective Hamiltonian H e ff. The sign of 
g also changes, but as noted above this is equivalent to 
replacing <v with —dJ. Thus this change has no effect on 
the eigenenergies, and the eigenstates only experience a 
change in the sign of those components where the plaque- 
tte is occupied by an odd number of dressed molecules. 
In addition, if particle-hole exchange changes the num- 
ber of atoms by N, then the eigenenergies are shifted by 
fA. 

From this investigation, we have seen that even on a 
single plaquette, the Hamiltonian -ff e // exhibits a num- 
ber of different types of ground state configurations, pos- 
sessing various forms of rotational symmetry (s-wave, d- 
wave, etc.). There are level crossings between these dif- 
ferent types of ground states as the detuning A is var- 
ied. The change of the ground state symmetry from s- 
wave to d-wave as one scans the parameter A from neg- 
ative to positive regions may be a general feature and 
not limited to a single plaquette. For a large lattice, this 
symmetry change might correspond to a quantum phase 
transition from the s-wave to the d-wave superfluidities 
|l9| . The states found in this work on a single plaquette 
also provide some basic entries for constructing the effec- 
tive many-body Hamiltonian for atoms in a quasi-two- 
dimensional optical lattice through the contractor renor- 
malization method 0, When the average filling 

number of the lattice is close to a half with hole doping, 
one expects that the basic degrees of freedom from each 
plaquette are the ground state configurations specified in 
Sec. II, the fermionic hole excitations given by the states 
in Sec. V, the bosonic hole-pair excitations specified in 
Sec. IV, and the bosonic spin excitations given by the 
states in Sec. Ill and II. The effective many-body Hamil- 
tonian will then describe the interaction between these 
basic degrees of freedom. So, it is our hope that the in- 
vestigation of the single-plaquette physics here will make 
it possible to better understand the physics of strongly 
interacting fermions on larger lattices. 

This work was supported by the NSF awards 
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